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ABSTRACT 

Building galaxy merger trees from a state-of-the-art cosmological hydrodynamics sim¬ 
ulation, Horizon- AGN, we perform a statistical study of how mergers and smooth accre¬ 
tion drive galaxy morphologic properties above z > 1 . More specifically, we investigate 
how stellar densities, effective radii and shape parameters derived from the inertia ten¬ 
sor depend on mergers of different mass ratios. We find strong evidence that smooth 
accretion tends to flatten small galaxies over cosmic time, leading to the formation of 
disks. On the other hand, mergers, and not only the major ones, exhibit a propensity 
to puff up and destroy stellar disks, confirming the origin of elliptical galaxies. We also 
And that elliptical galaxies are more susceptible to grow in size through mergers than 
disc galaxies with a size-mass evolution ro.s oc instead of ro.s oc 

depending on the merger mass ratio. The gas content drive the size-mass evolution due 
to merger with a faster size growth for gas-poor galaxies ro.s oc than for gas-rich 
galaxies ro.s oc Mg. 

Key words: galaxies: formation — galaxies: evolution — galaxies: interactions — 
galaxies: kinematics and dynamics — methods: numerical 


1 INTRODUCTION 

While our present understanding of galaxy evolution de¬ 
rives mainly from the nearby (z < 1) Universe, the bulk 
of today’s stellar mass formed around the broad peak of 
cosmic star formation history at z ^ 2 (e.g. Madau et al. 
1998; Hopkins & Beacom 2006). Although it represents a 
significant epoch in the evolution of the observable Uni¬ 
verse, the properties of galaxies remain largely unexplored at 
this epoch, as it has only recently become accessible by cur¬ 
rent observational facilities (CANDELS, GOODS, Herschel, 
ALMA, Chen et al. 2014; Lin et al. 2010; Cooper et al. 
2012; Lopez-Sanjuan et al. 2013). 

As a result, what constitutes arguably the most impor¬ 
tant aspect of hierarchical galaxy formation and evolution 
is still being debated. To what extent mergers, as opposed 
to secular evolution driven by (cold) gas inflows, explain 
the diversity of galaxies? The significance of mergers, con¬ 
sidered a cornerstone of the bottom-up growth of galaxies, 
has been heavily debated in recent work. They are certainly 

^ E-mail: welker@iap.fr 


capable of inducing star formation, black hole growth and 
morphological transformations (e.g. Springel et al. 2005), 
but it is not obvious that, at z > 1, mergers drive the 
evolution of galaxy properties like stellar mass, size and 
morphology (Shankar et al. 2004; Law 2009; Kaviraj et al. 
2013) considering the steady input of accreted streams of 
gas (Keres et al. 2005; Ocvirk et al. 2008; Dekel et al. 2009) 
and the gas-rich nature of galaxies (Tacconi et al. 2010; 
Santini et al. 2014). 

Observational studies suggest that large fractions of 
star-forming galaxies around z ^ 2 are not mergers 

but show kinematics and visual morphologies that are 
more consistent with systems dominated by turbulent 
discs (e.g. Forster Schreiber et al. 2006; Shapiro et al. 2008; 
Genzel & Burkert 2008; Mancini et al. 2011; Kaviraj et al. 
2013). In addition, many primordial spheroids that are form¬ 
ing the bulk of their stellar mass at z ~ 2 do not show 
the tidal features that would be expected from recent ma¬ 
jor mergers (Kaviraj et al. 2013). Notwithstanding these ad¬ 
vances, large statistical studies remain difficult, both due to 
the fact that observational efforts rely on pencil-beam sur¬ 
veys that are susceptible to low-number statistics and cosmic 
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variance but also that they are based on techniques that can 
differ significantly from study to study. Moreover, studies of 
galaxy merging at these redshifts are further complicated by 
the fact that normal star forming discs becomes more tur¬ 
bulent and asymmetric at earlier times, making them diffi¬ 
cult to separate from genuine mergers (Kaviraj et ah 2014; 
Huertas-Company et ah 2014). 

On the other hand, semi-analytical models and numeri¬ 
cal simulations propose that mergers can account for the size 
increase of local early-type galaxies if they are mostly dry 
(gas poor) and minor mergers (Boylan-Kolchin et al. 2006; 
Khochfar & Silk 2006; Mailer et al. 2006; Naab et al. 
2006, 2007; Bournaud et al. 2007; De Lucia 2007; 

Guo & White 2008; Hopkins et al. 2009; Nipoti et al. 
2009; Feldmann et al. 2010; Shankar et al. 2013; 
Bedorf & Portegies Zwart 2013). Dry minor mergers explain 
the loss of compactness of massive ellipticals at z < 2, where 
they are thought to take over smooth accretion processes 
in terms of stellar mass increase rates (Oser et al. 2010; 
Lackner et al. 2012; Hirschmann et al. 2012; Dubois et al. 
2013; Lee & Yi 2013). The dryness of low-redshift galax¬ 
ies is ensured either by the environment (for satellites 
infalling in groups and clusters) or by the presence of a 
supermassive black hole (BH) at the center of massive 
galaxies which powers feedback from the active galactic 
nuclei (AGN) (Di Matteo et al. 2008; Booth & Schaye 
2009; Dubois et al. 2012a). Together, these mechanisms 
allow for the formation of extended elliptical galaxies that 
would otherwise remain compact discs (Dubois et al. 2013; 
Ghoi et al. 2014). 

Multiple numerical studies also focused on a few ide¬ 
alised high resolution merger events to determine their 
impact on the morphology of the stellar component of 
galaxies (Bournaud et al. 2004, 2005; Naab & Trujillo 2006; 
Peirani et al. 2010). They found that while major merg¬ 
ers, or multiple minor mergers of stellar disks tend to pro¬ 
duce elliptical-like remnants, either disky or boxy depend¬ 
ing on the amount of gas available (Gretton et al. 2001; 
Naab & Burkert 2003; Naab et al. 2006; Qu et al. 2011), 
single minor mergers did not systematically destroy the 
primary disk but only thickened it (Quinn et al. 1993; 
Walker et al. 1996; Velazquez V White 1999; Younger et al. 
2007). On the other hand, the steady input of cosmo¬ 
logical gas accretion is able to rebuild the disc of galax¬ 
ies (Brooks et al. 2009; Agertz et al. 2009; Pichon et al. 
2011). Hence one needs to asses the relative impor¬ 
tance of mergers versus smooth accretion driven by 
the cosmic environment and to study its induced mor¬ 
phological diversity. With the advent of large-scale al¬ 
beit fairly well resolved cosmological hydrodynamical 
simulations such as Horizon-AGN (Dubois et al. 2014, 
see as well Devriendt et al. 2010; Khandai et al. 2014; 
Vogelsberger et al. 2014; Schaye et al. 2015 for similar sim¬ 
ulations performed with different numerical techniques), it 
has recently become feasible to investigate these different 
physical processes in detail and with sufficient statistics, a 
necessary requirement to truly unravel the impact of galaxy 
environment on their properties. 

Following up on Welker et al. (2014) which explored the 
fundamental role of mergers and smooth accretion on the 
evolution of the galactic spin and its orientation with re¬ 
spect to the cosmic web, we investigate the comparative role 


of mergers and smooth accretion of both gas and stars on 
defining and modifying the size and morphology of galaxies 
(see Fig 1). In the first section we describe the main fea¬ 
tures of the simulation and present the numerical analysis 
methods used to derive our results. In section 3, we evaluate 
the accretion rates of different types of mergers as well as 
smooth accretion over cosmic history. In section 4, we ana¬ 
lyze the impact of both processes on the growth of galaxies 
in the cosmic web, with specific emphasis on the different 
role played by gas and stars dominated mergers (equivalent 
to the dry/wet dichotomy used in lower z studies). In section 
5, we explore the competitive effects of smooth accretion and 
mergers on the morphology of galaxies and their correlation 
to the disk and spheroid abundances over the duration of 
the peak of cosmic star formation. Our main results and 
conclusions are summed up in section 6. 


2 NUMERICAL METHODS AND 
DEFINITIONS 

We briefly describe the cosmological hydrodynamical simu¬ 
lation used for analysis in this paper, Horizon-AGN, which 
is already described in more details in Dubois et al. (2014). 
The Horizon-AGN simulation is run with a AGDM cosmol¬ 
ogy with total matter density Qm = 0.272, dark energy den¬ 
sity = 0.728, amplitude of the matter power spectrum 
as = 0.81, baryon density fib = 0.045, Hubble constant 
Ho = 70.4kms“^ Mpc“^, and Us = 0.967 compatible with 
the WMAP-7 data (Komatsu 2011). The size of the simula¬ 
tion box is Lbox = 100 h~^ Mpc on a side, and the volume 
contains 1024^ dark matter (DM) particles, corresponding 
to a DM mass resolution of Mum, res = 8 x lO'^M©. The 
simulation is run with the RAMSES code (Teyssier 2002), 
and the initially coarse 1024^ grid is adaptively refined 
down to Ax = 1 proper kpc, with refinement triggered 
in a quasi-Lagrangian manner: if the number of DM par¬ 
ticles becomes greater than 8, or the total baryonic mass 
reaches 8 times the initial DM mass resolution in a cell. It 
lead to a typical number of 6.5 x 10® gas resolution ele¬ 
ments (leaf cells) in the Horizon-AGN simulation at z — 1. 
Heating of the gas from a uniform UV background takes 
place after redshift Zreion = 10 following Haardt & Madau 
(1996). Gas can cool down to 10^ K through H and He colli¬ 
sions with a contribution from metals using rates tabulated 
by Sutherland & Dopita (1993). Star formation occurs in re¬ 
gions of gas number density above no = 0.1 H cm“^ following 
a Schmidt law: p* = where p* is the star formation 

rate mass density, pg the gas mass density, e* = 0.02 the 
constant star formation efficiency, and tff the local free-fall 
time of the gas. Feedback from stellar winds, supernovae 
type la and type H are included into the simulation with 
mass, energy and metal release. The simulation also follow 
the formation of black holes (BHs), which can grow by gas 
accretion at a Bondi-capped-at-Eddington rate and coalesce 
when they form a tight enough binary. BHs release energy in 
a quasar/radio (heating/jet) mode when the accretion rate 
is respectively above and below one per cent of Eddington, 
with efficiencies tuned to match the BH-galaxy scaling rela¬ 
tions at z = 0 (see Dubois et al. 2012b, for details). 

Galaxies are identified with the most massive sub¬ 
node method (Tweed et al. 2009) of the AdaptaHOP halo 
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Figure 1. Rest-frame color images (w, g and i filters) of a sample of Horizon-AGN disk galaxies caught during their pre-merger phase at 
z = 2.2 (two leftmost columns) and their post-merger phase at 2 ; = 1.9 (two rightmost columns). The first and third columns are edge-on 
views, and the second and fourth columns are face-on views. Extinction by dust is not taken into account. Each frame is 100 kpc on a 
side. Ms is the stellar mass of the main progenitor in Mq units and 6m is the mass ratio of the merger (see text for exact definition). 
This figure illustrates the ability of mergers (major but also minor) to turn disk-like galaxies into spheroids. 


finder (Aubert et al. 2004) operating on the distribution of 
star particles with the same parameters than in Dubois et al. 
(2014). Unless specified otherwise, only structures with a 
minimum of Nmin = 100 particles are considered, which typi¬ 
cally selects objects with masses larger than 2x 10 ® M©. Cat- 
alogues containing up to ~ 150 000 galaxies are produced for 
each redshift output analysed in this paper (1.2 < z < 5.2). 
Note that, although sub-structures may remain, these popu¬ 
lations of galaxies are largely dominated by main structures. 

The galaxy catalogues are then used as an input to 
build merger trees with TreeMaker (Tweed et al. 2009). Any 
galaxy at redshift Zn is connected to its progenitors at red- 
shift Zn-i and its child at redshift Zn+i- We build merger 
trees for 22 outputs from z = 1.2 to z = 5.2 equally spaced in 
redshift. On average, the redshift difference between outputs 
corresponds to a time difference of 200 Myr (range between 
100 and 300 Myr). We reconstruct the merger history of 
each galaxy (halo) starting from the lowest redshift z and 
identifying the most massive progenitor at each time step as 


the galaxy or main progenitor^ and the other progenitors as 
satellites. Moreover, we check that the mass of any child con¬ 
tains at least half the mass of its main progenitor to prevent 
misidentifications. Note that such definition of mergers (vs 
smooth accretion) depends on the threshold used to identify 
objects as any object below the chosen threshold is discarded 
and considered as smooth accretion. 

We sort mergers in three categories depending on 
the mass fraction Sm = mniergers(zn-i^n)/Ms(zn), where 
mmergers(zn-i^n) is the stcllar mass accreted through merg¬ 
ers between Zn-i and Zn and Ms(zn) the stellar mass of the 
merger product at Zn- Major mergers are defined as mergers 
with Sm > 20%, minor mergers as mergers with 9% < Sm < 
20% and very minor mergers with 4.5% < Sm < 9%. Any 
merger with Sm < 4.5% is discarded and counted as smooth 
accretion. These bins are defined so as to be consistent with 
observational definitions of mergers using pairs of interact¬ 
ing galaxies, for which the observed mass ratio R is defined 
as = Msateiiite/Afgaiaxy and whcrc the subscripts indicate 
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secondary and main progenitors respectively, as defined in 
the previous paragraph of this section. Our bins thus corre¬ 
spond to = 1 : 4, = 1 : 10 and = 1 : 20. 

In order to preserve completeness, we define a second 
threshold and exclude galaxies with Mg < 5 x 10® M© from 
the galaxy sample used in our analysis. Satellites, however, 
are allowed to be less massive, in order for us to still cap¬ 
ture the whole range of merger ratios for the smallest of 
our galaxies. The galaxy threshold is identihed as the green 
vertical line on Fig. 2, which displays the mass function for 
all the structures identihed in the simulation (for compari¬ 
son we have also plotted the mass function of the Horizon- 
no AGN simulation, which is the same simulation performed 
without BHs and, therefore, AGN feedback). As can be seen 
on the hgure our sample is complete down to our strict se¬ 
lection threshold of Ms — 2x 10® M©, corresponding to 
galaxies with 100 star particles. Thus, the smallest mergers 
detectable for galaxies at the Mg < 5 x 10® M© threshold 
correspond to a mass ratio of 8m — 4.5%, which means that 
our merger classihcation is complete for our galaxy sample. 

The gas content and its properties (density, metallicity, 
pressure, temperature) of each galaxy is extracted from the 
AMR grid, considering all cells within its effective radius, 
(see Appendix B for the technical subtleties). 

Since the gas needs to be cold and dense enough to be 
eligible to form stars, we dehne “cold” gas (in the sense of 
star forming gas) the cells with a gas density higher than n > 
0.1 Hcm“® and a temperature T ^ 10^ K. We also define the 
gas fraction /gas of a galaxy as /gas = Mgas'^/(Mo .5 + 
with Mgas'^ the mass of cold gas and Mo. 5 the mass of stars, 
both enclosed within the sphere of radius ro. 5 . As can be 
seen in Fig. B2, this quantity decreases with stellar mass and 
redshift due to star formation and feedback, older galaxies 
becoming more massive after they used the gas available 
to form stars and/or after it has been blown out of them 
by AGN/supernova feedback. This evolution is consistent 
with previous numerical studies (e.g. Dubois et al. 2012b; 
Popping et al. 2014). 


3 MERGERS AND SMOOTH ACCRETION 
ACROSS THE PEAK OF COSMIC STAR 
FORMATION HISTORY 

While it is now well established that mergers have a signih- 
cant impact on z < 1 early-type galaxy sizes and kinematics 
(see references in introduction), it is not yet clear whether 

(i) this extends to the galaxy population at high redshift and 

(ii) over which timescale they are of importance, as many 
galaxies may not merge at all for long periods of time. Ob¬ 
servations of local galaxies not only suggest that early-type 
galaxies increased their size by 3-5 from z ^ 2, but also that, 
while most massive ones (Ms > 1.5x10^^ M©) roughly dou¬ 
bled their size from z ^ 1, smaller ones underwent a more 
limited growth between z ^ 1 and z =0 (by a factor 1.1 to 
1.3, Huertas-Gompany et al. 2013). From these results, we 
expect a growth by at least a factor 2 — 2.5 between z ~ 2 —3 
and z = 1 (Nipoti et al. 2012). To quantify the relative con¬ 
tribution of mergers and smooth accretion to the total mass 
budget of galaxies over the range of redshifts correspond¬ 
ing to the peak of cosmic star formation history down to 
z = 1 , we therefore compute the rates of galaxies having 



Figure 2. Galaxy stellar mass function in Horizon-AGN, for 
2 ; = 4 to ^ = 1.3. N is the number density of galaxies, Mg the 
stellar mass (together with Horizon-noAGN for comparison). The 
sharp cut-off at Mg = 10® M© corresponds to our completeness 
detection threshold. Observational points from CANDELS-UDS 
and GOOD-S surveys are rescaled from best fits in Mortlock et al. 
(2015) and Grazian et al. (2014) and overplotted. While mass 
functions are consistent at the high mass end, Horizon-AGN over¬ 
shoots the low-mass end by about a factor 3 in this redshift range. 
The vertical green line shows the selection threshold for our main 
progenitors candidates, chosen to enable us to completely track 
their mergers with galaxies up to 20 times smaller. 

undergone at least a merger within our mass fraction bins 
at these redshifts. 

We find that, at z = 1.2, around 35% of galaxies with 
Ms > 10^® M© have undergone at least one major merger, 
80% a minor merger, and 85% a very minor merger. These 
results are consistent with findings by Kaviraj et al. (2014), 
(our minor merger rates are slightly inferior due to a coarser 
redshift sampling). Fig. 3 (left panel) presents the evolu¬ 
tion of those rates with redshift, focusing on the sub-sample 
of galaxies in this mass range who possess a progenitor at 
z = 5.2 (sub-sample of 15 000 galaxies). It displays the evo¬ 
lution of the fraction of this sub-sample which remains free 
from mergers of a given type (major, minor and very mi¬ 
nor) as a function of cosmic time. It shows that over this 
4 Gyr period, ~ 50% of the sample undergoes a major 
merger and therefore that mergers, especially minor ones, 
are quite frequent over the whole redshift range. The sam¬ 
ple is affected by mergers at an average rate around 1 — 
2 X 10“® Gyr-^ h® Mpc“® and 3 - 5 x 10“^ Gyr”^ h® Mpc“® 
for all mergers and major mergers respectively. Note that 
these values are completely consistent with observations by 
Lotz et al. (2011) and in good agreement with the cumula¬ 
tive merger rates per galaxy derived from the Illustris sim¬ 
ulation (Rodriguez-Gomez et al. 2015) . 

The right panel of Fig. 3 focuses on galaxies which have 
had at least a merger between z = 5.2 and z = 1.2. It shows 
the probability distribution function (PDF) P(n, > Sm) for 
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Figure 3. Left panel: evolution of a sub-sample of galaxies identified at 2 ; = 1.2 with Ms > 10^° Mq and which can be tracked to z = 5.2. 
This panel shows the probability for galaxies not to undergo a major or minor merger during the redshift interval. Right panel: PDF 
of the number nm of galaxy mergers of a given mass ratio, (5m, undergone between redshifts 5.2 ^ z ^ 1.2. This PDF is restricted to 
galaxies with at least one very minor merger. Vertical lines show the average value for each sample. It illustrates the paucity of major 
mergers: most galaxies which merged, have had at most one major merger across this cosmic time interval, while they go through on 
average 2 to 3 mergers. Note the large variations from one galaxy to another. 


these galaxies to have undergone a number n of mergers of 
a given mass ratio, (5m, between z = 5.2 and z = 1.2. This 
PDF indicates that, while most of these galaxies underwent 
at most a single major merger, on average they undergo two 
to three merger events. 

Fig. 4 reveals that both mergers and smooth accretion 
ought to be taken into account to attempt to understand the 
morphology distribution of galaxies. This figure shows the 
evolution of the smooth accretion and merger contributions 
to the mass budget of galaxies across the peak of cosmic star 
formation history. The first panel presents the evolution of 
Am/Ms averaged for all galaxies with Ms > 10 ^° M© at 
2 : = 1 . 2 , with Am the mass accreted between two successive 
time steps (i.e. over a period of ~ 200 Myr) and Ms the 
stellar mass. The red and green curves correspond to the 
mass fraction accreted via smooth accretion of gas (green 
curve) and stars (red curve: gas+stars; i.e. including merg¬ 
ers with Ms < 2 X 10^ M© galaxies), and the blue curve 
corresponds to the mass fraction accreted through mergers. 
While at high-redshift {z ~ 5) young and small galaxies un¬ 
dergo a rapid relative mass growth through accretion of gas 
and swift merging of very small structures close to the detec¬ 
tion threshold, this activity settles around z ~ 3 — 4, when 
effects of smooth accretion and mergers on mass growth be¬ 
come comparable, until mergers slightly take over around 
z ~ 1.5. The net result is that at z = 1.2, ~ 45% of the 
galactic stellar mass can be attributed to in situ formation 
from smooth accretion of gas, as can be seen on Fig. 4 (sec¬ 
ond panel) . 

In conclusion, mergers and smooth accretion contribute 
equivalently to the galactic mass budget over the peak of cos¬ 
mic star formation history. It therefore seems that in order 
to understand the evolution of galactic sizes and morpholo¬ 


gies over this period, one needs to account for the possibility 
that these two processes play different roles. This is what we 
explore in the next sections. 


4 GROWTH OF GALAXY SIZES 

Fig. 5 shows the evolution of the stellar density, obtained 
by adding the masses of ah star particles enclosed within 
the half mass radius of the galaxy. Since the shape of 
galaxies can vary signihcantly over the cosmic time inter¬ 
val spanned by our study, especially when galaxies merge, 
we take anisotropy into account. More specihcahy, the den¬ 
sity p is dehned as p = 3 Mo. 5 /( 47 ra 5 c) with a > b > c 
the lengths of the semi-principal axes of the galaxy derived 
from the eigenvalues of the inertia tensor, as explained in 
Appendix A, and Mo. 5 is the sum of ah the masses of the 
star particles contained within its half mass radius. The left 
panel of the hgure shows the PDF of the relative density 
growth p = 2(y9n+i - pn)/{pn f- pn+i) , whcrc pn is the av¬ 
erage density of the galaxy within its half mass radius at 
time step n, stacked for each time output of the simulation 
between 1.2 ^ ^ 5.2. Notice how mergers tend to widen 

the distribution, populating the high-compactions and high- 
dilatation tails of the distribution. Looking at this panel, one 
might think that smooth accretion and very minor mergers 
tend to lower the density of the merger remnant on average, 
while minor and major mergers tend to increase it, but it is 
actually highly dependent on the initial morphology of the 
galaxy. This can be seen on the middle and right panels of 
Fig. 5. Galaxies that are initially disks show an increased 
stellar density after mergers, the effect being stronger the 
higher the mass ratio of the merger. On the other hand. 
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Figure 4. Left panel: evolution of Am/Mg over cosmic time for galaxies with Mg > 10^^ M© at 2 ; = 1.2, where Am is the mass increase 
due to in situ formed stars (green), merger with a companion (blue), or in situ formed stars combined with the ’diffuse’ accretion of 
stars (i.e. stars not identified as belonging to any galaxy, red) between two time steps (~ 200 Myr), and Mg is the stellar mass. We plot 
average values for all selected galaxies at every time output. Note that mergers and smooth accretion contribute similarly to the mass 
growth of galaxies from 2 ; = 3 onwards. Right panel: evolution of Mg^ooth/^s over cosmic time from 2 ; = 5.2 to 2 ; = 1.2, where Mg^ooth 
is the mass of stars either produced in situ from the gas component, or accreted ’smoothly’ (i.e. star particles not associated with a 
galaxy above our mass resolution threshold), for galaxies with Mg > 10^® Mq at 2; = 1.2. At 2; = 1.2, about half of the stellar mass of 
these galaxies comes from such smooth accretion processes. 


galaxies which begin as spheroids tend to have their stellar 
density decreased by mergers, the effect being statistically 
stronger for minor mergers. 55% of spheroids which merge 
betray a decrease in stellar density after ~ 200 Myr, and 
only 40% of the non-merger galaxies do the same. It is inter¬ 
esting to notice that, in both cases, this increase/decrease 
in average density is related to how skewed the distribu¬ 
tion becomes and not only to a global drift towards posi¬ 
tive/negative values. As a result, minor and major mergers 
of spheroids are much more likely to trigger important de¬ 
creases in stellar density (by more than a factor 2) than 
smooth accretion: 16% of cases versus 5% respectively, and 
even as low as 3% if events where stars are accreted below 
our galaxy mass threshold are discarded. Around 8% of ma¬ 
jor mergers and even fewer minor mergers trigger dilatations 
by more than a factor 5. 

Similarly, 73% of disk galaxies increase their density af¬ 
ter merging (against 63% for non-mergers). However, and 
more importantly, 30% of these mergers increase it by at 
least a factor 2 as compared to only 9% for smooth accre¬ 
tion. Finally, only 10% of the major mergers and the minor 
mergers lead to compactions by more than a factor 5. 

These results statistically support the claim that merg¬ 
ers turn disks into denser structures while they tend to lower 
the density of spheroids. Moreover, although major mergers 
are found to be quite rare, the ability of the more frequent 
minor mergers to trigger effects of comparable amplitude 
points towards an important role of multi-minor mergers 
in driving the size-mass relationship of early-type galax¬ 
ies (Kaviraj 2014). 

We further analyze the role of mergers in driving galaxy 


stellar density evolution by looking at the relation between 
growth in stellar mass and growth in stellar half-mass ra¬ 
dius. To do so, we compute the evolution of the logarithmic 
derivative of the half mass radius ro .5 with respect to Mg as 
a function of the mass ratio 5m — Am/Mg oc Alog^^o^s? 
where Am is the stellar mass accreted between two con¬ 
secutive outputs through mergers (blue curve) or smooth 
accretion (red and green curves) for all galaxies and all time 
outputs between z = 5.2 and z = 1 . 2 . Note that Mg or Mo .5 
are equivalent for the purpose of this comparison so we use 
them interchangeably. The result is shown in the top left 
panel of Fig. 6 . While mergers and smooth accretion drive a 
similar amount of mass growth (see the left panel of Fig. 4 
in the previous section), mergers are much more efficient 
drivers of galaxy size growth. Ignoring the very shallow de¬ 
pendence on 5m^ smooth accretion processes lead to an av¬ 
erage radius-mass relation ro .5 oc Mf with a = 0.1 ±0.05. 
Slightly higher values of a are reached for higher 5m (see 
left panel of Fig. 6 ), but always remain within a factor 2, 
i.e. a ^ 0.2. This dependence on 5m can be explained by the 
fact that higher mass ratio values most often correspond to 
lower stellar masses, for which cold flows bring in more spe¬ 
cific angular momentum than hot phase accretion in more 
massive galaxies (Kimm et al. 2011). For mergers, we ob¬ 
tain ro .5 oc Mf with [3 = 0.85 ± 0.3. /3 also increases with 
5m up to values ~ 1.2. These values are consistent with 
observations (e.g. Newman et al. 2012; Cimatti et al. 2012; 
Huertas-Company et al. 2013; van der Wei et al. 2014, who 
found a value of /3 around ~ 0.6 — 0.8 for early-type galax¬ 
ies), and together with the smaller values of a, support 
the idea that the size growth of galaxies is mostly driven 
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Figure 5. Left panel: PDF of the density growth ratio /i = 2(pn+i — pn)l{pn +Pn+i) , where pn is the density of the galaxy within its 
half mass radius at time step n, for different merger mass ratios. This ratio is calculated for galaxies with Mg > 10®-^Mq over each time 
step between 1.2 ^ 2 ; ^ 5.2, and all these timesteps are then stacked. Each vertical dashed line shows the average value for the merger 
mass ratio bin of the corresponding color. Mergers have a tendency to widen the distribution and increase the stellar density, especially 
major mergers (vertical dashed line on the positive side of p values). However, this behavior is actually different for galaxies which are 
disks prior to the merger (for which the stellar density rises: middle panel) and for those which are originally spheroids (for which the 
stellar density decreases: right panel). 


by mergers (Boylan-Kolchin et al. 2006; Nipoti et al. 2009; 
Feldmann et al. 2010; Dubois et al. 2013). The bottom pan¬ 
els show the evolution of the average half-mass radius ro .5 
as a function of redshift when we split our galaxy sample in 
bins of ’’final” mass (i.e. galaxy masses at redshift z = 1 . 2 ; 
bottom left panel) and in bins of constant stellar mass (i.e. 
independent of redshift; bottom right panel). Once again, 
the results shown on Fig. 6 are consistent with the over¬ 
all evolution of the size-mass relationship from observations 
such as Huertas-Company et al. (2013) (though our simu¬ 
lated galaxies are a factor of ~ 2 larger): galaxies of a given 
stellar mass display much larger radii at 2 : = 1.2 than their 
counterparts of similar stellar mass at z = 5.2, most of the 
growth taking place between z = 3 and z = 1.2. More specif¬ 
ically, at 2 : = 1.2 galaxies with Mg > 10^°'^ M© display 
an average half-mass radius twice to 3 times bigger than 
their counterparts at 2 : = 3. (see bottom right panel), while 
one can see on the bottom left panel that galaxies reaching 
Mg > 10^°'^ M© at 2 ; = 1.2 have also seen their half-mass 
radius grow by a factor 2 to 3 since 2 : = 3, and by a factor 
4 since 2 : = 5.2. 

However, Fig. 6 does not distinguish gas-rich and 
gas-poor mergers. This could be potentially important 
as gas-poor mergers are known to trigger intense size 
growth of local early type galaxies (e.g. Naab et al. 2007; 
Feldmann et al. 2010), whereas accretion of gas (by gas-rich 
mergers or smooth accretion) is thought not to be able to 
since their gas shocks radiatively and loses angular momen¬ 
tum, therefore piling up in the central region of the galaxy 
where it rapidly turns into stars and causes size contrac¬ 
tion. The right panel of Fig. 7, lends statistical support to 
this claim. It shows the average value of the merger mass 
ratio 6m which leads to a given relative variation of the half 
mass radius Aro.s/ro.s, for mergers with /gas > 0.6 (blue 
data points) and /gas < 0.6 (yellow data points). From this 


data, one can see that radius contraction (negative values 
of Aro.s/'To.s) is confined to gas rich minor mergers (blue 
data points with 0.09 < Sm < 0.2). The corresponding yel¬ 
low data points for gas poor mergers are below - or very 
close to - the smooth accretion threshold (lower horizon¬ 
tal dashed line), indicating that smooth accretion of gas 
is in fact the leading advection process in those cases. In¬ 
terestingly enough, major mergers Sm > 0.2 statistically 
never lead to a compactification of galaxies, regardless of 
whether they are gas rich or not: the violent disruption that 
they occasion does not translate into a funelling of mate¬ 
rial to the central region as it does for minor mergers, but 
as an extended redistribution of it. Note that the threshold 
of /gas = 0.6 is chosen high compared to the values tradi¬ 
tionally used to define wet and dry mergers at low redshift 
because galaxies are more gas rich on average in the redshift 
range of this study. One can get an idea of how much smaller 
the sample gets when this threshold is lowered to /gas = 0.2 
by looking at the left panel of Fig. 7. 

This panel presents the dependence of the logarith¬ 
mic derivative of the half mass radius ro.5 on the mass 
ratio Sm for our sample of galaxies split into different 
pre-merger gas fraction bins. One can see that star rich 
mergers with /gas < 0.6 (in yellow), especially minor ones 
(0.09 < Sm < 0 . 2 ) induce a more efficient radius growth 
than their gas dominated counterparts of similar mass ratio 
(in blue). Gas deprived mergers with /gas < 0.2 (pink curve), 
whether major or minor, lead to a rapid growth of the effec¬ 
tive radius compatible with ro .5 cx: M/ where 7 = 2 zb 0.5. 
This power law index is in excellent agreement with predic¬ 
tions from Hilz et al. (2012), and consistent with previous 
numerical studies (Boylan-Kolchin et al. 2006; Nipoti et al. 
2009; Feldmann et al. 2010), although slightly higher, which 
therefore lends extra support to the scenario involving mul- 
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Figure 6. Top left panel: evolution of the logarithmic derivative of the half-mass radius ro.5 with respect to the mass ratio 5m = Am/Ms 
(Am is the stellar mass gained between two consecutive outputs) of the merger or smoothly accreted material for all galaxies and all 
time outputs between 2 ; = 5.2 and 2 ; = 1.2. Filled blue symbols indicate mergers, red and light green ones represent smoothly accreted 
mass, including a stellar component (red) or gas only (green). Whilst the evolution is linear in each case, the dependence of radius growth 
on the mass ratio is found to be much steeper for mergers. Top right panel: same plot as the top left panel, except we have split the 
merger sample according to different pre-merger morphologies: disks (green) and spheroids (yellow). The steepness of the radius versus 
5m relation appears mainly caused by minor merger disruption of the disks. Bottom left panel: evolution of the half-mass radius as a 
function of redshift for galaxies split into bins of different mass at redshift 2 : = 1.2. Bottom right panel: same plot as in the bottom left 
panel but for galaxies split into redshift independent mass bins. 


tiple dry mergers to explain the loss of compacity of massive 
early-type galaxies at low redshifts. 

Going back to the top right panel of Fig. 6, we see that 
the dependence of the size-mass relationship on merger mass 
ratio can be interpreted as a morphological effect: galaxies 
that are spheroids prior to the merger (yellow data points) 
systematically grow in size almost indistinctively with mass 


ratio (except for the most extreme major mergers), whilst 
disks (green data points) exhibit a size growth proportional 
to the accreted mass ratio over the same range in 5m. Note 
that van der Wei et al. (2014) find a different size-mass evo¬ 
lution for early-type and late-type galaxies with [3 ~ 0.75 
and 13 ~ 0.22 respectively with negligible evolution with 
redshift. In our simulation, we find that spheroids (i.e. early- 
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Figure 7. Left panel: evolution of the relative variation of the half-mass radius as a function of the merger mass ratio for all galaxies 
which undergo a merger between 2 ; = 5.2 and 2 ; = 1.2 and for different pre-merger gas fractions: /gas < 0.2 (pink symbols), /gas < 0.6 
(yellow symbols) and /gas > 0.6 (blue symbols). The error bars plotted correspond to Icr errors. Horizontal dashed lines represent 
ro .5 oc Ms ^ with 7 = 1 and 7 = 2 which are predicted size-mass relations for (dry) major and minor mergers using the virial theorem 
(Hilz et al. 2012; Dubois et al. 2013) . Note how the presence of gas limits the radius growth. Right panel: average mass ratio versus 
relative variation of the half-mass radius Aro.s/ro.s for mergers with gas fraction /gas < 0.6 (yellow) and /gas > 0.6 (blue). Horizontal 
dashed lines show major/minor/smooth accretion separation thresholds in 5m. The vertical dashed line indicates the border between 
expansion (positive values) and contraction (negative values). Note how radius contraction is confined to wet minor mergers. 


type galaxies) have /3 = 1.2 on average and disks (i.e. late- 
type galaxies) have /3 ~ —0.5 for low values of 5m and 
/3 ~ 0.5 for large values oi 5m^ which shows a similar dis¬ 
crepancy of the size-mass evolution between different galaxy 
morphologies to that observed in van der Wei et al. (2014). 
This stresses the need to study the morphology of our galax¬ 
ies in further detail, and we now turn to this issue. 


5 IMPACT ON GALAXY MORPHOLOGIES 

Focussing on more accurate morphological parameters, 
Fig. 8 displays the time evolution - for galaxies which do 
not merge - of the cumulative PDFs of the principal semi¬ 
axis ratios = c/a, ^2 — c/b and — b/a with a > b > c 
of the inner half-mass stellar component, derived by calcu¬ 
lating the inertia tensor of the galaxy (see Appendix A). 
One can see from this figure that, while ^3 tends to remain 
constant over cosmic time, with a value strongly peaked at 
1 (large axis equals to intermediate axis), both average val¬ 
ues of and ^2 decrease at an average rate of almost 10% 
per Gyr, from 0.64 and 0.74 down to 0.54 and 0.64 respec¬ 
tively in the 4 Gyr which separate z = 5.2 and z = 1.2. 
For reference, note that an infinitely thin and homogeneous 
disk has = <^2 = 0 and ^3 = 1 while a perfect sphere has 
= ^2 = ^3 = 1- Our result indicates that smooth accretion 
and consecutive in situ star formation tend to flatten galax¬ 
ies over time along the minor axis, which coincides with the 
spin axis. However, this accretion has no significant effect 
on the circularity of galaxy disks. 

This morphological transformation strongly depends on 


the mass and morphology of galaxies. As explained in Ap¬ 
pendix A, we define spheroids as galaxies with > 0.7 and 
^2 > 0.8 and disks as galaxies with < 0.45 and ^2 < 0.55. 
Fig. 9 displays the evolution of the principal semi-axis ra¬ 
tio PDFs for galaxies classified as spheroids (blue curves 
and symbols) and disks (red curves and symbols) pre-merger 
(galaxies are excluded from the sample when they merge), 
and for two different mass bins. The upper panels focus on 
galaxies with a stellar mass comprised between 10^'^Mq and 
lOio.5^0, lower panel on more massive galaxies with 
Ms > 10^°'^ M©. This mass threshold corresponds to the 
transition mass above which galaxies embedded in filaments 
decouple from their vorticity quadrant and display a spin 
perpendicular to their closest filament (see Dubois et al. 
2014), and also to the transition in gas accreted onto the 
galaxy between cold and hot mode (e.g. Dekel & Birnboim 
2006; Ocvirk et al. 2008). The figure reveals that the de¬ 
crease rate in and ^2 is much faster, around 20% per 
Gyr for spheroids with masses below the transition mass (or 
from average values of = 0.7 and ^2 = 0.8 to 0.56 and 
0.66 respectively). On the other hand, disks tend to thicken 
slightly on average (going from = 0.46 and ^2 = 0.56 to 
0.5 and 0.64 respectively). This behaviour for the disks at 
least partially arises from the limited maximum spatial res¬ 
olution of the simulation (1 kpc). By definition disks with 
scale heights below this value are artificially ’puffed up’ to 
1 kpc and any accretion of new material, no matter how 
dynamically cold, can only result in increasing this minimal 
numerical scale height. This especially alters the shape of 
small galaxies, for which the scale length is also poorly re¬ 
solved. For galaxies with masses above the transition mass 
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Figure 8. Cumulative PDFs of the principal semi-axis ratios for all galaxies with Mg > 10 ^'^Mq which do not merge. Colors represent 
evolution with cosmic time from 2 : = 5.2 (dark red) to 2 : = 1.2 (light orange) with an average time step of ~ 200 Myr. The Icr poissonian 
error bars are overplotted for all bins that have a non-zero probability. Smooth accretion tends to flatten galaxies over time. 


(bottom panels of Fig 9), we do not observe any significant 
impact of smooth processes on the morphology indicators 
which remain constant on average. 

Interpreting our results in the light of the scenario 
described by Codis et al. (2012) and Laigle et al. (2015), 
whereby small galaxies acquire their spin through angu¬ 
lar momentum transfer from the vorticity quadrant they 
are embedded in, this ’’flattening” effect can be under¬ 
stood as the (re)-formation of disks in high vorticity re¬ 
gions at the heart of cosmic web hlaments. In other words, 
smooth accretion tends to (re)-ahgn galaxies with their 
nearest hlament (Tillson et al. 2012; Welker et al. 2014; 
Danovich et al. 2014; Pichon et al. 2014) where the domi¬ 
nant component in this process, for galaxies below the tran¬ 
sition mass, is coherent gas feeding from cold flows. At the 
opposite end of the mass spectrum, galaxies above the tran¬ 
sition mass accrete material from multiple quadrants and/or 
smaller amounts of material along a unique hlament. The 
angular momentum streamed to the core of these massive 
galaxies from multiple directions is more likely to cancel 
out, which results in little to no effect of smooth accretion 
on the morphology of the galaxy. These results reinforce ear¬ 
lier hndings that the underlying cosmic web plays a major 
role in shaping galaxy properties. 

As can be seen on Fig. 10 (and D2), mergers trig¬ 
ger very different evolutions for disk galaxies (respectively 
spheroids). This hgure showcases a qualitative difference: 
contrary to smooth accretion, both minor and major merg¬ 
ers strongly change the galaxy morphology, leading to much 
more spheroidal/ehiptical structures. However, this only oc¬ 
curs for disks: as can be seen in Fig. D2, for galaxies ini¬ 
tially identihed as spheroids, minor mergers behave more 
like smooth accretion, flattening the galaxy whereas major 
mergers preserve, by and large, their morphology. Looking 
at Fig. 10, there are quantitative differences between minor 
and major mergers. Whilst major mergers clearly destroy 
disks (the average value of the PDF of — cja shifts from 
0.45 to 0.62, red and dark blue curves respectively), minor 
mergers have a more limited effect (^1 PDF average value 
shifted from 0.45 to 0.52 only). Finally, the effect of very 
minor mergers (light orange curve) is closer to a thickening 


of the disks than an actual destruction of them and an al¬ 
teration of the galaxy morphology. It is important to notice 
that ah mergers also trigger an increase in the scatter of the 
distribution of galactic disks morphology indicators as 
the slope of the PDFs becomes shallower for mergers than 
smooth accretion. This effect is stronger for higher merger 
mass ratios. 

These hndings corroborate the view that major and 
multiple minor mergers can lead to galaxies with simi¬ 
lar morphologies, destroying disks and turning them into 
spheroids (e.g. Bournaud et al. 2007). This allows to over¬ 
come the tension occasioned by the paucity of major merg¬ 
ers: minor mergers are much more frequent events, allow¬ 
ing for the formation of a much larger spheroid population. 
An illustration of this phenomenon is given in Fig. 1 which 
depicts rest-frame false color images of a sample of disk 
galaxies in the Horizon-AGN simulation before and after 
major/minor merging observed through the u, g and i Al¬ 
ters. 


6 CONCLUSION 

Galaxy growth in the cosmic web involves a wide range of 
processes from anisotropic accretion to supernovae and AGN 
feedback whose effect can either add up or cancel one an¬ 
other, resulting in the observed diversity in morphologies, 
kinematics and colors of galaxies. While the interplay be¬ 
tween these phenomena is undoubtedly complex, the ap¬ 
proach we implemented in this work which consists in fo¬ 
cussing on a couple of well dehned processes (smooth accre¬ 
tion against mergers) and identify their impact - whether 
re-enforcing or competing - on specihc galactic properties 
(size, morphology) still yields some interesting results: 

• Mergers and smooth accretion augment galaxy masses 
across the peak of cosmic star formation history, in amounts 
that are statistically comparable. As a result, at 2 ; = 1.2, 
galaxies with Mg > 10^° M© have acquired 55% of their stel¬ 
lar mass via smooth accretion and 45% via mergers. How¬ 
ever, while smooth accretion is a steady process with regular 
impact on stellar mass over cosmic history, mergers are vi- 
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Figure 9. Cumulative PDFs of the principal semi-axis ratios for galaxies which do not merge. Colors evolve with cosmic time from 
dark red (z = 5.2) to light orange {z = 1.2) for disks and navy blue (z = 5.2) to light blue (z = 1.2) for spheroids, with an average 
time step of ~ 200 Myr. The Icr poissonian error bars are overplotted for all bins that have a non-zero probability. Top panel: galaxies 
with 9.5 < log(Ms/MQ) < 10.5. Bottom panel: galaxies with log(Ms/M 0 ) > 10.5 (spin flip and cold/hot mode for accretion transition 
mass: see text for detail). Smooth accretion has little impact on the morphology of galaxies above the transition mass but clearly flattens 
galaxies below it. 



Figure 10. Cumulative PDFs of the morphology indicators for different merger mass ratios and for galaxies with Mg > 10^® Mq 
identifled as pre-merger disks. The poissonian Icr error bars are indicated for all bins with more than 10 galaxies. The results are stacked 
for each time output between 2 ; = 5.2 to 2 ; = 1 . 2 . Mergers broaden the morphology distribution and quantitatively destroy disks. This 
effect strengthens with increasing merger mass ratios. 
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olent processes which occur on average twice in the history 
of a galaxy over this epoch. 

• Mergers and smooth accretion augment galaxy sizes 
across the peak of cosmic star formation history, especially 
major mergers, but this growth strongly depends on red- 
shift and gas fraction. We also found that while mass is 
accreted, the mean density also rises for galaxies which are 
pre-merger disks, suggesting a gravitational contraction dur¬ 
ing the merger phase, while the inverse is true for pre-merger 
spheroids which on average expand after merging. 

• For mergers of mass ratio 5m^ the relative increase in 

radius is found to evolve as a power law of the stellar mass 
ro.5 = j^o.65+(5m smoothly accreted material of com¬ 

parable mass ratio proves to be much less efficient in growing 
galaxy radii ro.5 = . Moreover, while the growth of 

spheroid sizes shows little dependence on 5m (ro.5 oc Mg), 
- even for smallest minor mergers which is consistent with 
the idea that material is then smoothly accreted within the 
galactic plane -, disks show a stronger dependence on 5m^ 
even contracting trend when subjected to minor mergers. We 
interpret this result as the destruction of disks and redistri¬ 
bution of their stellar component in a more tightly packed 
spheroidal volume, which causes the effective half mass ra¬ 
dius to decrease even though the amount of mass accreted 
actually increases. 

• Gas fraction also plays an important part in determin¬ 
ing the size growth consecutive to mass accretion. As ex¬ 
pected, gas dominated mergers induce a much more lim¬ 
ited growth in size than star dominated ones. In such gas 
rich mergers, the remnant appears to be more compact. We 
interpret this as the result of gas shocking, losing angular 
momentum and being transported to the central parts of 
the galaxy where it forms stars, seemingly triggering a grav¬ 
itational contraction of the galaxy. On the opposite, star 
dominated mergers (with /gas < 0.2) induce an increased 
growth in radius with no significant dependence on merger 
mass fraction, but the steepest dependence on stellar mass 
that we measure (ro.5 = ^s)- 

• These accretion processes are found to have a strong im¬ 
pact on galaxy morphologies. Smooth accretion tends to flat¬ 
ten small galaxies along their spin axis, consistent with the 
idea that those galaxies are embedded in a vorticity quad¬ 
rant of cosmic filament which feeds them angular momentum 
coherently along the filament direction. This effect is even 
clearer for the sub-sample of spheroids fed by this smooth 
accretion which evolve to resemble the disk population in 
just over 2 Gyr. In contrast, mergers tend to destroy disks 
and form spheroids, except for very minor mergers - which 
only thicken them -, in agreement with the idea that in this 
case the satellite is slowly stripped from its gas and stars 
in the galactic plane of the main progenitor. But our main 
result is that minor mergers are responsible for a compara¬ 
ble amount of disk destruction than major mergers, coupled 
with a strong contraction effect when the minor merger hap¬ 
pens to a gas-dominated (/gas >0.6) galaxy. 

• These results altogether statistically favour a scenario 
whereby galaxies grow their stellar mass by smooth accre¬ 
tion of gas, in situ formation and mergers in comparable 
amounts, but grow in size mostly through merging: disk 
(gas-dominated) galaxies merge to become more compact 
spheroids while spheroids lose their compactness through 
these same minor mergers. Occasionally, dramatic growth in 


size through rare major mergers and multiple, gas-deprived 
minor mergers happens. Non-merging spheroids with masses 
up to a transition mass around M© then rebuild disks 

from coherent smooth accretion. Above this mass the coher¬ 
ence of streams is lost and morphology is preserved. 

Though this first study supports - in a full cosmological 
context using the Horizon- AGN simulation - the consistent 
galaxy growth model that has emerged from previous nu¬ 
merical studies of different types of mergers, further inves¬ 
tigation is required to extend these results down to z = 0 
and to specify in detail the role played by galactic physics - 
more specifically supernovae and AGN feedback - in shap¬ 
ing these results. Analyzing more specific merger parameters 
such as the impact parameter and the orbital-to-intrinsic 
angular momentum transfer rate will also be necessary to 
understand the scattering of the morphology and size dis¬ 
tributions induced by mergers and understand their overall 
impact on observed galaxies in the local Universe. Finally, 
the internal kinematics of the galaxy population also need 
to be examined more closely. 


ACKNOWLEDGMENTS 

This work was granted access to the HPC re¬ 
sources of CINES (Jade) under the allocation 
2013047012 and c2014047012 made by GENCL This 
research is part of the Spin(e) (ANR-13-BS05-0005, 
http: //cosmicorigin. org) and Horizon-UK projects. 
Let us thank D. Munro for freely distributing his Yorick 
programming language and opengl interface (available at 
http: //yorick. sourceforge. net/). We warmly thank 
S. Rouberol for running the Horizon cluster on which the 
simulation was post-processed. Part of the analysis was 
also performed using the DiRAC facility jointly funded 
by BIS and STEC. The research of JD is supported by 
Adrian Beecroft, The Oxford Martin School and STEC. CP 
thanks Churchill college for hospitality while this work was 
finalized. 

References 

Agertz O., Teyssier R., Moore B., 2009, MNRAS, 397, L64 
Aubert D., Pichon G., Golombi S., 2004, MNRAS, 352, 376 
Bedorf J., Portegies Zwart S., 2013, MNRAS, 431, 767 
Booth G. M., Schaye J., 2009, MNRAS, 398, 53 
Bournaud F., Gombes F., Jog G. J., 2004, A&A, 418, L27 
Bournaud F., Jog G., Gombes F., 2007, A&A, 476, 1179 
Bournaud F., Jog G. J., Gombes F., 2005, A&A, 437, 69 
Boylan-Kolchin M., Ma G.-P., Quataert E., 2006, MNRAS, 
369, 1081 

Brooks A. M., Governato F., Quinn T., Brook G. B., Wad- 
sley J., 2009, ApJ, 694, 396 

Ghen G., Smail I., Swinbank A., et al. 2014, ArXiv e-prints 
Ghoi E., Ostriker J. P., Naab T., Oser L., Moster B. P., 
2014, ArXiv e-prints 

Gimatti A., Nipoti G., Gassata P., 2012, MNRAS, 422, L62 
Godis S., Pichon G., Devriendt J., Slyz A., Pogosyan D., 
Dubois Y., Sousbie T., 2012, MNRAS, 427, 3320 
Gooper M. G., Griffith R. L., Newman J. A., Goil A. L., 
Davis M., Dutton A. A., Eaber S. M., Guhathakurta P., 


The rise and fall of z > 1 disk galaxies 13 


Koo D. C., Lotz J. M., Weiner B. J., Willmer C. N. A., 
Yan R., 2012, MNRAS, 419, 3018 
Cretton N., Naab T., Rix H.-W., Burkert A., 2001, ApJ, 
554, 291 

Danovich M., Dekel A., Hahn O., Ceverino D., Primack J., 
2014, ArXiv e-prints 

De Lucia G., 2007, Cosmic Frontiers ASP Conference Se¬ 
ries, 379, 257 

Dekel A., Birnboim Y., 2006, MNRAS, 368, 2 
Dekel A., Birnboim Y., Engel G., Freundlich J., Goerdt 
T., Mumcuoglu M., Neistein E., Pichon C., Teyssier R., 
Zinger E., 2009, Nature, 457, 451 
Devriendt J., Rimes C., Pichon C., Teyssier R., Le Borgne 
D., Aubert D., Audit E., Colombi S., Courty S., Dubois 
Y., Prunet S., Rasera Y., Slyz A., Tweed D., 2010, MN¬ 
RAS, 403, L84 

Di Matteo T., Colberg J., Springel V., Hernquist L., Siiacki 
D., 2008, ApJ, 676, 33 

Dubois Y., Devriendt J., Slyz A., Teyssier R., 2012a, MN¬ 
RAS, p. 2188 

Dubois Y., Devriendt J., Slyz A., Teyssier R., 2012b, MN¬ 
RAS, 420, 2662 

Dubois Y., Gavazzi R., Peirani S., Silk J., 2013, MNRAS, 
433, 3297 

Dubois Y., Pichon C., Welker C., et ah 2014, MNRAS, 444, 
1453 

Feldmann R., Carollo C. M., Mayer L., Renzini A., Lake 

G. , Quinn T., Stinson G. S., Yepes G., 2010, ApJ, 709, 
218 

Forster Schreiber N. M., Genzel R., Lehnert M. D. e., 2006, 
ApJ, 645, 1062 

Genzel R., Burkert A. e., 2008, ApJ, 687, 59 

Grazian A., Fontana A., Santini P., Dunlop J. S., Ferguson 

H. C., et ah 2014, ArXiv e-prints 

Guo Q., White S., 2008, MNRAS, 384, 2 
Haardt F., Madau P., 1996, ApJ, 461, 20 
Hilz M., Naab T., Ostriker J. P., Thomas J., Burkert A., 
Jesseit R., 2012, MNRAS, 425, 3119 
Hirschmann M., Naab T., Somerville R. S., Burkert A., 
Oser L., 2012, MNRAS, 419, 3200 
Hopkins A. M., Beacom J. F., 2006, ApJ, 651, 142 
Hopkins P. F., Bundy K., Murray N., Quataert E., Lauer 
T. R., Ma C.-R, 2009, MNRAS, 398, 898 
Huertas-Company M., Kaviraj S., Mei S., O’Connell R. W., 
Windhorst R., Cohen S. H., Hathi . P., Koekemoer A. M., 
Licitra R., Raichoor A., Rutkowski M. J., 2014, ArXiv 
e-prints 

Huertas-Company M., Mei S., Shankar F., Delaye L., Rai¬ 
choor A., Covone G., Finoguenov A., Kneib J. P., Le F. O., 
Povic M., 2013, MNRAS, 428, 1715 
Kaviraj S., 2014, MNRAS, 440, 2944 
Kaviraj S., Cohen S., Ellis R. S. e., 2013, MNRAS, 428, 
925 

Kaviraj S., Cohen S., Windhorst R. A., Silk J., O’Connell 
R. W., Dopita M. A., Dekel A., Hathi N. P., Straughn A., 
Rutkowski M., 2013, MNRAS, 429, L40 
Kaviraj S., Devriendt J., Dubois Y., Slyz A., Welker C., 
Pichon C., Peirani S., Le Borgne D., 2014, ArXiv e-prints 
Kaviraj S., Huertas-Company M., Cohen S., Peirani S. e., 
2014, MNRAS, 443, 1861 

Keres D., Katz N., Weinberg D. H., Dave R., 2005, MN¬ 
RAS, 363, 2 


Khandai N., Di Matteo T., Croft R., Wilkins S. M., Feng 
Y., Tucker E., DeCraf C., Liu M.-S., 2014, ArXiv e-prints 
Khochfar S., Silk J., 2006, AJ, 648, L21 
Kimm T., Devriendt J., Slyz A., Pichon C., Kassin S. A., 
Dubois Y., 2011, ArXiv e-prints 
Komatsu E. et al ., 2011, ApJ Sup., 192, 18 
Lackner C. N., Gen R., Ostriker J. P., Joung M. R., 2012, 
MNRAS, 425, 641 

Laigle C., Pichon C., Codis S., Dubois Y., Le Borgne 
D., Pogosyan D., Devriendt J., Peirani S., Prunet S., 
Rouberol S., Slyz A., Sousbie T., 2015, MNRAS, 446, 2744 
Lambas D. G., Maddox S. J., Loveday J., 1992, MNRAS, 
258, 404 

Law D. R., 2009, PhD thesis, California Institute of Tech¬ 
nology 

Lee J., Yi S. K., 2013, ApJ, 766, 38 

Lin L., Cooper M. C., Jian H.-Y., Koo D. C., Patton D. R., 
Yan R., Willmer C. N. A., Coil A. L., Chinch T., Croton 
D. J., Gerke B. F., Lotz J., Guhathakurta P., Newman 
J. A., 2010, ApJ, 718, 1158 

Lopez-Sanjuan C., Le Fevre O., Tasca L. A. M., Epinat 
B., Amram P., Contini T., Garilli B., Kissler-Patig M., 
Moult aka J., Paioro L., Perret V., Queyrel J., Tresse L., 
Vergani D., Divoy C., 2013, A&A, 553, A78 
Lotz J. M., Jonsson P., Cox T. J., Croton D., Primack J. R., 
Somerville R. S., Stewart K., 2011, ApJ, 742, 103 
Madau P., Pozzetti L., Dickinson M., 1998, ApJ, 498, 106 
Mailer A., Katz N., Keres D., Dav R., Weinberg D., 2006, 
ApJ, 647, 763 

Mancini C., Forster Schreiber N. M., Renzini A., Cresci 
G. e, 2011, ApJ, 743, 86 

Mortlock A., Conselice C. J., Hartley W. G., Duncan K., 
Lani C., et ah 2015, MNRAS, 447, 2 
Naab T., Burkert A., 2003, ApJ, 597, 893 
Naab T., Jesseit R., Burkert A., 2006, MNRAS, 372, 839 
Naab T., Johansson P., Jeremiah P.and Efstathiou G., 
2007, ApJ, 658, 710 

Naab T., Khochfar S., Burkert A., 2006, ApJ, 636, L81 
Naab T., Trujillo L, 2006, MNRAS, 369, 625 
Newman A. B., Ellis R. S., Bundy K., Treu T., 2012, ApJ, 
746, 162 

Nipoti C., Treu T., Auger M. W., Bolton A. S., 2009, ApJ 
Let., 706, L86 

Nipoti C., Treu T., Leauthaud A., Bundy K., Newman 
A. B., Auger M. W., 2012, MNRAS, 422, 1714 
Ocvirk P., Pichon C., Teyssier R., 2008, MNRAS, 390, 1326 
Oser L., Ostriker J. P., Naab T., Johansson P. H., Burkert 
A., 2010, ApJ, 725, 2312 

Peirani S., Crockett R. M., Geen S., Khochfar S., Kaviraj 
S., Silk J., 2010, MNRAS, 405, 2327 
Pichon C., Codis S., Pogosyan D., Dubois Y., Desjacques 
V., Devriendt J., 2014, ArXiv e-prints 
Pichon C., Pogosyan D., Kimm T., Slyz A., Devriendt J., 
Dubois Y., 2011, MNRAS, pp 1739—h 
Popping G., Somerville R. S., Trager S. C., 2014, MNRAS, 
442, 2398 

Qu Y., Di Matteo P., Lehnert M. D., van Driel W., 2011, 
A&A, 530, AlO 

Quinn P. J., Hernquist L., Fullagar D. P., 1993, ApJ, 403, 
74 

Rodriguez-Gomez V., Genel S., Vogelsberger M., Sijacki D., 
Pillepich A., Sales L. V., Torrey P., Snyder G., Nelson D., 


14 C. Welker et al. 


Springel V., Ma C.-P., Hernquist L., 2015, ArXiv e-prints 

Santini P., Maiolino R., Magnelli B., et al. 2014, A&A^ 562, 
A30 

Schaye J., Crain R. A., Bower R., et al. 2015, MNRAS, 
446, 521 

Shankar F., Marulli F., Bernardi M., Mei S., Meert A., 
Vikram V., 2013, MNRAS, 428, 109 

Shankar F., Salucci P., Granato G. L., De Zotti G., Danese 
L., 2004, MNRAS, 354, 1020 

Shapiro K. L., Genzel R., Forster Schreiber N. M. e., 2008, 
ApJ, 682, 231 

Springel V., Di Matteo T., Hernquist L., 2005, ApJ Let., 
620, L79 

Sutherland R. S., Dopita M. A., 1993, ApJ Sup., 88, 253 

Tacconi L. J., Genzel R., Neri R., et al. 2010, Nature, 463, 
781 

Teyssier R., 2002, A&A, 385, 337 

Tillson H., Devriendt J., Slyz A., Miller L., Pichon C., 2012, 
ArXiv e-prints 

Tweed D., Devriendt J., Blaizot J., Colombi S., Slyz A., 
2009, A&A, 506, 647 

van der Wei A., Franx M., van Dokkum P. G., et al. 2014, 
ApJ, 788, 28 

Velazquez H., White S. D. M., 1999, MNRAS, 304, 254 

Vogelsberger M., Genel S., Springel V., Torrey P., Sijacki 
D., Xu D., Snyder G., Nelson D., Hernquist L., 2014, MN¬ 
RAS, 444, 1518 

Walker I. R., Mihos J. C., Hernquist L., 1996, ApJ, 460, 
121 

Welker C., Devriendt J., Dubois Y., Pichon C., Peirani S., 
2014, MNRAS Letters, 445, L46 

Younger J. D., Cox T. J., Seth A. C., Hernquist L., 2007, 
ApJ, 670, 269 


APPENDIX A: INERTIA TENSOR 

The inertia tensor lij of a galaxy is computed from its star 
particle distribution (indexed by 1) and calculated at the 
centre of mass of the galaxy, according to the definition: 
lij = T>im\Sij— x\.x^j), where mi is the mass of 
star particle I and x{ its position in the barycentric coordi¬ 
nate system of the galaxy. As a 3x3 real symmetric matrix, 
the inertia tensor can be diagonalized, with its eigenvalues 
Al > A2 > A3 being the moments of inertia relative to the 
basis of principal axes ei, 62 and 63. The lengths of the 
semi-principal axes a, b and c (with a > b > c) are straight¬ 
forwardly derived from the moments of inertia: 

a — (5/Mo.5)VAi + A 2 - A3 , along e3, 

5= (5 /Mo. 5 )\/Ai + A 3 - A2 , along 62, 
c = (5 /Mo. 5)VA3 + A2 - Al, along ei . 

Therefore the morphology of a galaxy can be easily 
interpreted using its axis ratios: ^1 = c/a, ^2 = c/b and 
^3 = b/a. A perfectly round and infinitely thin disk has 
^1 = 0, ^2 = 0 and ^3 = 1. For a Milky Way like galaxy 
(including stars from the inner bulge+thin disk+bar), one 
gets ^1 = 0.06, ^2 = 0.07 and ^3 = 0.98. The limited spatial 
resolution (1 kpc) of the Horizon-AGN simulation, prevents 
us from obtaining disks as thin as these. We therefore iden¬ 
tify disks with our most flattened ellipsoids. More specif- 



Figure Bl. 2D sketch of the gas cell assignment procedure for 
one galaxy from Horizon-AGN (shown as a face-on projected gas 
density map). The thick red circle represents the effective radius 
ro.5 around the galactic center of mass and the white squares 
the AMR grid with different levels of refinement. The green tick 
indicates when a cell or sub-cell is counted as belonging to the 
galaxy (see text for detail). 


ically, we adopt < 0.45 and ^2 < 0.55 as a definition 
for disks and > 0.7 and ^2 > 0.8 to define spheroids. 
Other galaxies are simply classified as ellipsoids. Though 
these cuts may appear a crude approximation, they are ac¬ 
tually quite consistent with 3-D axis ratios reconstructed 
from observations (Lambas et al. 1992). Note that we also 
define the morphology of our galaxies using star particles en¬ 
closed within their half mass radius sphere. We found that 
this is more robust than using all the star particles identified 
by our halo finder especially for post-merger remnants, as 
these can exhibit elongated tidal features which persist for 
a considerable amount of time. 


APPENDIX B: GAS CONTENT 

Fig. Bl displays a sketch of the systematic procedure used 
to extract gas content from AMR cells for all galaxies in 
Horizon-AGN. 

Fig. B2 shows the evolution of the average gas fraction 
across the peak of cosmic star formation history for galaxies 
of different masses. Our results are consistent with previ¬ 
ous numerical investigations (e.g. Popping et al. 2014). /gas 
decreases with redshift as star formation consumes and feed¬ 
back expels the available gas. 

For each galaxy in the sample, we define the maximum 
radius rmax as the distance between the galactic center of 
mass (GOM) and the furthest star particle, the effective ra¬ 
dius ro.5 as the half stellar mass radius and Arceii = ro.s/lO. 
AMR cells with a size dceii larger than Arceii are subdivided 
in 2^^^^ sub-cells with ric such that dsub-ceii < Arceii- AMR 
cells counted as belonging to the galaxy are: 1) AMR cells 
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Figure B2. Evolution of the gas fraction /gas in the redshift 
range 1.8 < 2 ; < 4.5 for different mass bins, where M 0.5 is the 
stellar mass enclosed within the half mass sphere and Mgun the 
mass of the sun. Results are consistent with previous simulations 
(e.g. Popping et al. 2014). /gas decreases with redshift as star 
formation consumes the available gas and/or feedback blows it 
out of the galaxies. 

with a size dceii < Arceii and a center within the sphere of 
radius ro.s centred on the COM 2) sub-cells of larger AMR 
cells with a length dceii < Arceii and a center within the 
sphere of radius ro.s centred on the galaxy COM. This pro¬ 
cedure is illustrated on Fig. B1 which shows a 2D sketch of 
the cell selection process on a face on projected gas density 
map for a galaxy from Horizon-AGN with a post-merger 
sub-structure at z = 3. 


APPENDIX C: MASS 

Fig. Cl shows the evolution of the average stellar mass of 
galaxies with Mg > 10^° M© at z = 1.2, across the peak 
of cosmic star formation history. Note that the stellar mass 
growth of galaxies remains steady for most of the evolution. 
The knee at redshift z = 1.5 corresponds to a peak in the 
merger rate and smooth accretion observed at the same red¬ 
shift (see Fig. 4), which is due to the extra level of refinement 
added at this particular redshift. 



Figure Cl. Evolution of the average stellar mass of galaxies 
with Ms > 10 ^*^ M 0 at 2 ; = 1 . 2 . 



Figure Dl. Cumulative PDF of — cja over one time step, with 
all timesteps between 2 ; = 5.2 and 2 ; = 1.2 stacked, for different 
merger mass ratios and for galaxies with Mg > IO^^Mq. 


APPENDIX D: MORPHOLOGY 

Fig. Dl shows the cumulative probability distributions of 
the morphologic ratio = c/a over one time step, with all 
timesteps between 2; = 5.2 and 2; = 1.2 stacked, for different 
merger mass ratios. This is plotted for all galaxies regardless 
of their morphology. Fig. D2 displays the same quantity for 
galaxies identified as spheroids, i.e. with morphologic ratios 
> 0.7 and ^2 > 0.8, before they merge. 
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Figure D2. Same as Fig. D1 but for galaxies which are classified 
as spheroids before they merge (see text for detail). 











